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1.0 SUMMARY 



In this essay a detailed analytical and experimental study of 
a new, high speed, linear torque motor has been carried out. First, the 
need for such a motor is shown. This is followed by a general description 
of the actual motor. The equations of motion and impedance for the motor 
are developed theoretically and analytical expressions for all the constants 
appearing in these equations are derived. Experimental values of all the 
constants are then found and experimental runs are made of the phase and 
amplitude characteristics, the impedance, and the transient response to a 
step input. A comparison is made of the theoretical and experimental 
results including a discussion of linear and non-linear regions. The 
effect of the motor in the closed servo loop is then discussed as well as 
the effect of current and voltage feedback. Finally, the conclusions 
arrived at as a result of the study are stated. 



2.0 INTRODUCTION 



2.1 STATEMENT OF PROBLEM 

At present there exists considerable interest in and need for 
small size, light weight, high speed, linear, torque motors at the Applied 
Physics Laboratory, The Johns Hopkins University, Silver Spring, Maryland. 
These motors are used in closed loop, electrohydraulic servomechanisms 
(see Fig. 1). To properly analyze these servomechanisms, the operating 
characteristics of all components should be known thoroughly. It is the 
object of this thesis to conduct a detailed analytic and experimental 
study of one of these high speed torque motors. 

The motor was developed independently by the Applied Physics 
Laboratory and Curtiss-Wright for servosystem usage. A similar type 
was constructed by the Dynamics and Control Laboratory of the Massachu- 
setts Institute of Technology. Two types of motors were built for the 
Applied Physics Laboratory — one by Curtiss-Wright Corporation and the 
other by Bendix Corporation (Eclipse-Pioneer Division). No thorough 
combined analytical and experimental studies have been conducted up to 
present due to the lack of time. 

This analysis is conducted on the motor built by Eclipse- 
Pioneer since it is used in preference to the Curtiss-Wright motor in 
which the phase shift varies greatly with frequency because of large 
eddy current losses in the magnetic structure. 

2.2 GENERAL DESCRIPTION OF MOTOR 

The motor is essentially a polarized relay consisting of two 
magnetic yokes energized by a permanent magnet, a moving armature made 
of laminated iron, and two armature coils. The armature is pivoted at 
Point P, midwa^ between the laminated pole faces of the yoke, as shown 
in Fig. 2. One end of the armature extends out beyond the gaps and is 
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connected to the load. 



The permanent magnet produces a constant flux, 0^, which divides 
into two paths and generates equal flux densities at each of the four gaps 
when the armature is centered. If, under these conditions, currents are 
allowed to flow through the two armature coils then each current produces 
a flux. The difference of these two fluxes will be called the armature 
flux and will be considered to be produced by the difference of the two 
currents actually flowing or what can be called the differential current. 
This armature flux causes the armature to become polarized in such a way 
as to be attracted by one pair of diagonal poles and repelled by the 
other. There is a torsional spring restraining the armature motion which 
together with the design of the magnetic circuit causes the physical dis- 
placement of the armature to be proportional to the current in the coils. 
Reversing the direction of the differential current in the coils reverses 
the direction of motion. Torque is produced by simultaneously raising the 
current in one coil and lowering it in the other. The motor is operated 
with a quiescent current flowing in each coil. The motor operates with 
direct current power and is in principle similar to a Class A amplifier. 

An exploded view of the motor is shown in Fig. 3 and the flux 
paths for the armature centered and off-center in Figs. UA and UB . 



FIELD 

LAMINATIONS 




FIG. 3. EXPLODED VIEW OF LINEAR TORQUE MOTOR BUILT BY BENDIX, ECU PSE- PION EER 
FOR THE APPLIED PHYSICS LABORATORY, THE JOHNS HOPKINS UNIVERSITY. 
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3.0 THEORY 



In this section the theoretical equations of motion and impedance 
will be derived and theoretical values or expressions for all the constants 
appearing in them. 

3.1 DERIVATION OF EQUATIONS OF MOTION AND IMPEDANCE 

3.11 Equivalent Circuits and Equations of Motion and Impedance 
3. HI Linear relationship between force produced by the motor 
and the differential armature current. 

The first thing to be found is the relationship between the force 
exerted by the armature arm and the differential armature current. The 
flux in the armature, 0 , is proportional to the difference of the currents, 

a 

i, flowing in the armature coils, i.e. 0 = a i, where a is a constant as 
long as the iron is not saturated. In Fig. 2 it can be seen from the geome- 
try of the relay that the gap lengths g-^-f- g 2 always equal g^ 4- g^. The 

permanent flux, 0 , emanating from the permanent magnet is considered con- 
P 

stant. It will be assumed that as long as the movement of the armature 
is small and the permeance of the gaps low, the flux through the armature 
from the permanent magnet is negligible. 

The force at any gap is given by the following formula according 
to Ref. 1, p. 202. 

2 

F = B S where F = force in lbs. 

7£ B = flux density, kilomaxwells/ sq. in. 

S = area of gap, sq. in. 

But, B = 0/S where 0 O is flux at any given gap 

so F _ i r 0j 2 S 

IS] 72 

If 0 a — 0 and the armature is centered, there will be no force 
present to move it away from center, i.e. F-^= F 2 and F^ = F^ in Figo 5>. 



If 0, then the total flux at each gap will be (see Fig. 2) 

h = V 2 f 0 a /2 

0 ? = 0 p /2 - 0 a /2 
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The force at each gap will then be: 

F x = ( 0 p /2 4- 0 a /2) 2 S 

^2 72 



F 2 = ( 0/2 - 0 a /2)2 S 

—2 -2 ?2 



F 3 = ( 0 p /2 - 0 a /2) 2 S 
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/V (force exerted by valve) 



FIG. 5 FORCES ACTING ON ARMATURE. 

The forces acting on the armature are as indicated in Fig. 5 
(all attractive). 

Writing the torque equation for the armature, 

F,X + F, X * F„X + F_X f F (X + Y) 

1.4 2 3 v 



or F f F. — F - F 0 — F (X + I) 
1 U 2 3 v 



X 



Substituting for F^, ? 2 > F^, F^ the expressions just derived and at 
the same time putting in the numerical values of X and Y, namely X= 0.781 
in. and Y = 0.625 in,, one finds that 

F v= 

61i.7S 

or F 1 = b 0 where b = 0 _ 

V 3. p 

6U.7S 

But 0 =r a i 

'a 

therefore, F y = a b i = K^i where = ab 

will be known as the electromechanical coupling constant and will be 
discussed in more detail later on. It is now seen that if the assumptions 
made in this derivation hold, the force is proportional to the differential 
current. 

3.112 MECHANICAL SYSTEM 

Now, considering the motor as a mechanical system, any force 
produced by the current in the armature coil must be numerically equal 
to the forces due to the inertia of the armature plus the moving mass 
of the piston and piston rod, to viscous friction, and to an effective 
spring force. By effective spring force is meant the combination of the 
force due to the torsional centering spring and the force due to permanent 
flux in the armature as a result of displacement. This mechanical system 
can be represented diagrammatically as follows: 
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FIG. 6. DIAGRAMMATIC REPRESENTATION OF MECHANICAL SYSTEM 



The symbolic diagram of Fig 0 6 is: 




FIG. 7. SYMBOLIC DIAGRAM OF MECHANICAL SYSTEM 



Writing the differential equation for this system. 



F = M d 2 x f f dx f kx 
5 

dt dt 



F = 

v 

M = 



f = 
k = 

x - 



force produced, lbs. 
total effective moving mass, 

lb sec^/i- n - 

viscous friction, lb sec/in. 
effective spring coefficient, 
lb/in. 

armature displacement, in. 



Since F v = J^i, 

F v - K«i — M d^x -L f dx i- kx 

_ 

dt dt 



3.113 ELECTRICAL SYSTEM 



The electrical system for the motor may be represented as 

in Fig. 8. 
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e 

i 
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1 




R c = resistance of coil, ohms. 

L = inductance of coil, henry. 

Kg = back e,m.f., volt sec/in. 
e = voltage feeding armature coil, 
volts. 



FIG 8. ELECTRICAL SYSTEM FOR THE MOTOR 



The output voltage of the last stage of the amplifier must feed the re- 
sistance and inductance of the armature coil and also overcome a back 



electromotive force due to the velocity of the armature. An equivalent 
circuit is drawn in Fig. 9. 






FIG. 9. EQUIVALENT ELECTRICAL CIRCUIT FOR THE MOTOR 

The voltage equation, by Kirchhoff's law, is: 

yue g = Ri -i- Ldi 4 - K g dx where R = r p 4- R c 

dt dt 

3.11U EQUATIONS OF MOTION OF THE MOTOR 

With these equations 

Li — Md^x f dx +. kx 
dt? dt 

and / JLe g- Ri +• Ldi 4 - Kgdx 
dt dt 

the output displacement per unit voltage input to the motor and also the 
impedance may be determined. Assuming that the system is initially at 
rest and using the Laplace transformation technique, one may write the 
following equations: 

(Ms^ +- fs 4 - k)x - K^i = 0 

(Kgs)x -+• (Ls -+- R)i =/ <e g 

A bar (-) over a letter indicates the Laplace transform of that quantity. 
Solving for x and i by means of determinants, 




r p = plate resistance, ohms. 

amplification factor. 

e — grid volts. 

6 









' 

I 






13 . 



0 



g 



-% 

Ls -I - R 



x 



% e g 
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%M e g 



LMs- 5 t- (fL 4- MR)s 4- (kL 4 fR +■ KgK^s +- kR 



i =• 



Ms 2 f fsfk 
K b s 



o 

/ 4e c 



(Ms* 4- fs 4- k)/<eg 



Ms fs k -K« 
Kgs Ls 4- « 



LMs 3 1- (fL r MR)s 2 +■ (kL 4- fR f Kg%)s 4- kR 



The Laplace transform of the displacement per unit volt input is then 






Me LMs^ 4- (fL -I- MR)s 2 +■ (kL 4- fR 4- KgK^s +- kR 



ify/LM 



4- (f/M r R/L)s ? 4- (k/M 4- fR/LM 4- %%/LM) 4- kR/LM 



upon dividing 



numerator and denominator by LM. 

The Laplace transform of the impedance^ Z, is e/i where e =jxe g when the 
plate resistance,!^, is zero, or in other words, one has a zero impedance 
source. 

Z = / ne g _ L j s 3 4- (f/M 4- R/L)s 2 4- (k/M 4- fR/LM 4- KgK M /LM)s 4- kR/LM] 
i s 2 4- (f/M)s 4* k/M 



where R now equals only R c . 

The Laplace transform of the displacement per unit ampere is found by 



dividing x by i and this reduces to 

_x_= % 

I (Ms 2 +- fs t- k) 

3.115 MORE EXACT EQUATIONS OF MOTION 

It should be noted that the equations of motion derived in 
3.HU are simplifications of the true situation. No attention has been 
paid to the characteristic of the connection from the motor armature to 
the piston valve. A symbolic diagram taking this into account is shown 
in Fig. 10. 
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FIG. 10. SYMBOLIC DIAGRAM OF MOTOR TAKING INTO ACCOUNT ELASTANCE OF 

CONNECTION FROM MOTOR TO VALVE. 

k v — spring coefficient of valve, 
kf — spring coefficient of connection. 

— displacement of piston rod. 
f v — viscous friction of piston valve. 

M^. — mass of valve and rod. 

M a -- mass of armature. 
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Writing the equations of motion, 



= (R c +- r p )i -t L di 4- % dx 
dt dt 



and Ifyi = M & d 2 x -4- kx 4 - Q 
dt? 



where Q = (x x - x) = f 4- k y x 1 



f rom which x _ 



f v*l + (k v - k l ,x lj 
k l 



A s'ngle dot over a letter indicates the first derivative with respect to 
time of that quantity and two dots the second derivative. 

Taking the Laplace transforms, 

Ae g — (R c t- Tp -4- Ls)i f Kgsx 

%i = (M a s 2 -t- k)x 4 - (M v s 2 4 - f v s k^^ =. (K a s 2 f- k - k-^x 4-k J x 1 



i = 1 (M a s 2 +- k - k^)x 4- k-^x^ 

% 



and x -_1 (M v s^ +- f y s +- 1 c,. - k^)x 1 
k. 



Let R — R c +- r p c 



Theny/e = (R 4- Ls)i 4 - K^sx 



or substituting for i and x 



/<e g = (R 4- Ls) 1 



% 



,M a s 2 k ^ — kf) j - 1_ (M y s f- f v 8 4 k y - k^x^-t-k^x^ 



4 



J 



- K b s 1 ( M v s 2 4- f^ 4- k y - k 1 )x 1 
*1 



Dividing through by x-^ and then taking the reciprocal, 
*1 1 



(R4- Ls) 1 


(M a s 2 4- k - k x )) - 1 (M v s 2 4 f v s f ^ - k x ) 


4- 
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-K B s (J^s i-f v s 


k m 


l k l 







4 -ky-ki) 



Multiplying the numerator and denominator by fyk^, we have 



% *1 



y^eg Ls ) 



(M a s 2 +- k - k iJf-l (M v s 2 +- f v s 4- k v - 

t k l 

-K B K M (M v 3 2 ^f v s +- - k^s 



-«K K ! 

A’s^-t- B’s 1 * f C’ S 3 4- D’s 2 f E'afF’ 



k l j+- k l| 



where A 
R* 
C* 



MaM„L 



D = R 



w +- w 

H a f y E +- L[M v (k - k x ) I- H a (k v - k^T- K b K„M v 
« v (k- k l> + M a (k v “ k J| +• k V k - k l> + Wm 



E* - f v R(k- k x ) 4- L(k-k 1 )(k v - k x ) r %%(k v - k x )- k x 2 L 



F* - (k-k 1 )(k v - k 1 )R-k 1 2 R 



When k^ is very large, i.e. k^ > oo 

A* — 0 

B* =•• 0 

C" = -L(M 4- M ) 

v a 

D* =-R(M v -h M a ) - Lf y 
E"“-f v R ~L(k v -h k) - K b K m 
F* = -(k 4- k y )R 

Substituting these values into the above equations and letting 
M = M v 4- M a and fy = f , 



JL- , «h 

/Ue LMs 3 -t- (fL -I- MR)s 2 4- (kL +■ f R +- KgK^s 4- kR 
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which is the same equation for displacement per unit input voltage derived 
in 3.11h for the simpler system assuming infinite elastance for the con- 
nection from the armature arm to the piston valve, an equivalent effective 
mass equal to the sum of the mass of the armature and valve, and that the 
viscous friction acts at the armature rather than in the valve itselfo 

3.12 Derivations of Mathematical Expressions for Constants Appearing 
in above Equations and Calculation of Set of Specific Values for 
Motor under Study . 

In this section the various constants appearing in the equations 
of motion and impedance will be derived mathematically. Later these theo- 
retical constants will be compared to those found experimentally. 

3»121 Effective Mass M 

As has already been noted the mass of the armature must be 
accelerated. Since this motor is to be used as a part of a system and it 
is the object of this thesis to determine its operating characteristics in 
that system an adjustment must be made in this term to take the moving mass 
of the piston rod and valve piston into account. For the purpose of sim- 
plification the mass of the armature can be combined in an equivalent mass 
with the total moving mass of the valve piston and rod. 

M = M e + M v 

M g = mass equivalent of armature 

M y = total mass of valve piston and piston rod 

(measured and found to be 17 grams = O°037!i> pounds) 

M = total effective moving mass converted to position 
of valve 

To find the equivalent mass of the armature, the moment of 
inertia of the armature was calculated from its dimensions and shape 
using the fact that the moment of inertia of a body is equal to the 
integral of the elemental mass multiplied by the square of the distance 



from this elemental mass to the axis about which the moment of inertia is 
sought and the parallel axis theorem. The calculation for I a , the moment 
of inertia of the armature, will be found in appendix 8.1. Theoretically, 

I a =• li8„12 gm-cm 2 . Therefore the equivalent mass is 

M Q - I a = U8.l/(3.57) 2 = 3.78 gram or 0,00835 pound. 

7 ? 

r — 3.57 cm is distance from pivot point P to point of 
application of valve force. 

Now the total effective moving mass, M, can be found. M — M e +- My 

M = 3.78 +■ 17.0 - 20.8 gm or 1.188x10^ lb sec 2 /in. 

3.122 Effective Spring Coefficient k 

For static equilibrium the forces on the armature arm must be 
equal, i.e. F y = k g x — K^i +- k^x 

where k = actual spring rate of torsional 
s 

centering spring. 

kf — anti-spring rate due to the permanent 
flux through the armature when armature 
is off center. 

= electromechanical coupling constant 

The force, k^x, due to the permanent flux flowing through the armature 
when it is off center (see Fig. UB ) causes the armature to move farther 
off center. This force plus the force due to the flow of armature current 
must equal the force due to the centering spring for static equilibrium. 

The displacement sensitivity per ampere, x/i, is found as follows: 

k g x - K^i 4- kfX 

x/i ^/(kg - k f ) 

If (k s — kf) is set equal to k, k may be looked upon as an effective spring 



coefficient. 



In order to find the value of k necessary to get mechanical 



resonance at, for example, 210 cycles per second, use is made of some well 
known relationships for second order systems. (These relationships may be 
found in any standard text on servomechanisms, for example, Ref. 2.) 

co =s 






angular velocity at undamped 
natural resonant frequency, 
rad/sec. 

k = effective spring coefficient, 
lbs/in. 

M = mass of system, lb sec^/in. 
or k= Man = (1.188x10”^) ( 13 20) 2 = 207 lbs/in. 

This particular resonant frequency was chosen because of the effective 
spring rate found in the experimental section and it must be borne in 
mind that it is only the true resonant frequency for a system having no 
damping. 



3*123 Viscous Friction f 

Assume J = 0.1 as a reasonable estimate as to the amount of 
damping present in the system, where J is defined as the ratio of the 
actual damping to the critical damping. 




These relationships -f _ r _ f f = viscous friction coefficient 

give J “ 2KCin 2pT 

or f = J 2 {W 

= (0.1) (2) I.l88xl0“^x207 =0.031)4 lb sec/in 

A useful parameter for calculations was found to be f/M which for J = 0„1 

was f _ O.Q31h = 26U/sec 

K Iol88xl0 4 ’ 



TABLE I. VALUES OF f/M CORRESPONDING TO DIFFERENT VALUES OF J 



I 


f/M 


0.018 


1 * 7.5 


0.100 


261 * 


0.113 


300 


0.300 


800 



The value of J will be used as a parameter to show the variation 
of impedance and output displacement per volt as the amount of damping in 
the system is varied. 



3.12U Inductance L 

L — N0 a 76^0 x 6 $60 8.U henries 

i x 10 8 0.060 x 10 8 

N = number of turns on armature coil 
0 a = armature flux in maxwells 
i = differential current in amperes 

The values of 0 a and i used to determine the theoretical value of inductance 
were taken from Fig. 55 in section 8.2„ 

3.12^ Back Electromotive Force K 

B 

This method of determining the back electromotive force (e.m.f.) 
is due to Mr. R. B„ Higley and Mr. J. Wheeler of Curtiss-Wright Corporation 
(see Ref. 3, p. 17). The leakage flux, 0^, however, has been determined in 
a different and more rigorous manner. 
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x = valve displacement, in. 

y = armature displacement 
in air gap, in. 



FIG. 11. RELATIONSHIP BETWEEN VALVE AND ARMATURE DISPLACEMENTS 



It can be seen from Fig. 11 that x _ y 

T7Ho5 oTTHI 



or y = 0-.731 = o.fo 

l . ii66 

When the armature is centered all gap lengths are equal or 



= &2 — £>3 — 

When the armature is not centered, g ? = gj = g — 0.6x 



Si = eu = e + °* 6x 



The permeance at any gap with the armature centered is 



(P= M&. 

g 



A>= permeability 
S = area 
g = gap length 



with the armature off center, (P, = (P )t = 

1 u ' g+ 0.6x 

and P ? = P~\ = A-S 












: 



Let us draw the magnetic circuit for the flux of the permanent magnet 
(see Fig. 12). 




arwi^k are 
reluctance of air gaps 
1) 2 , 3) and li <> 

ia reluctance of 
armature . 



FIG. 12. MAGNETIC CIRCUIT FOR PERMANENT MAGNET FLUX 



An expression for the flux, 0^, flowing through the armature due 
to the permanent magnet can be found by solving the three magnetic loop 
equations. These equations are: 

tfl+fyWl' — ^i 02' ~ (^2 03 / =NI 

~ < R\ 0i / ‘l" (^i ■+* < ^3 +C ?L)02 / ~ (R L 03 / = O 

— ^2 01 ' 02 /_ *“ ( < ^ 2 + '^ U '*" ^0^2 =0 

The flux we are seeking, 0^, is equal to ^6, — <^b / 

Actually solving we find 

0 L _ . NI 

^?L (^3 +• ^ U ) (^1 1~ ^ 2 ) -I - ^3^U (^1 +■ ^ 1^2 (^3 +■ ) 

or, in terms of permeances, remembering (P^=- (P^ and (P^— (P , 

4) 2 ^ + 1 ) + k Ik + ?) 



Dividing numerator and denominator by 

0L = £) 



(«r f w) 



/ 






-h 



(P,(P 3 



Putting the fractions over a common denominator and multiplying numerator 
and denominator by (P^ (P^ 



t L== ni(^ 3 -^,) 



&±JL , 2 



But we may justifiably assume, since the reluctance of the path through the 
armature is approximately zero, that the permeance,*^, is approximately 
infinite or certainly very much larger than^ 4“ (P^ „ 

Therefore, =0 

~^r~ 



or 



P = NI ((P 3 ~(Pl) 



Substituting f or v and (P-. , one ha 3 



Z( T _ NI 
L ~1 



A? _ 

'g-0.6x g +- 0.6x 



NI /, 3 



2(0o6)x 

g 2 -0.36x 2 



NI// 0 S 1.2 

T 7 



g 2 -0.36x 2 



Ibis checks with the expression given by Higley and Wheeler 
solving in this manner requires only one assumption, n amely. 



in Ref. 3, but 
that (P^)(P^V(P X 



It should be noted that one can reduce 0^, which is desired small, by either 
keeping x very small or by increasing g, the gap length, with the armature 
centered. 

To get the back e.m.f. due to 0^, use is made of the fact that 
0^ is a function of displacement x, and x is a function of time. Therefore, 
0^ is a function of time. 

E = - Na0LlG~® volts where N is number of turns of coil 

dt 0L is flux in maxwells 

E — — Nd0] ax 1CT 8 
dx dt 



E= dx where IC= — N d0L 10~ 8 volt sec/inch 

" dt ® dx 



If we differentiate 0 L =: NI/ x.S 1.2 



with respect to x, we get 



_g 2 -.36x^ 

d0 L= NI^S 1.2 ^(g 2 -.36x 2 )- 1 -f- 2(.36)x 2 (g 2 -.36x 2 )~ 2 J 



_ NI//.S 1.2 

T' 



(g 2 j- . 36x 2 ) 
(g ? -.36x 2 ) 2 



Substituting this into our expression for Kg, we get 




1.2 



g 2 -(- .36 x 2 1lo~ 8 



g -.36 x 2 ) 2 



volt sec/inch 



To use this formula the NI of the permanent magnet must be determined. 

The method used to determine this was as follows. The flux flowing in 
the armature due to the permanent magnet with the armature centered was 
assumed negligible. Also there was assumed to be no leakage flux, toith 
these conditions plus the further assumption that the reluctance of the 
iron path is negligible compared to the air gaps, all the magnetic poten- 
tial difference of the permanent magnet will appear across air gaps 1 and 2 



(see Fig. 2). 

r = tf?0 P 




£ 






^ = magnetic potential difference, 
gilberts 

ft — reluctance, gilbert/maxwell 
0p — permanent flux, maxwells 
/t, = permeability, maxwell/NI-cm 
Bp — flux density in any gap due to 
permanent magnet, gauss 



Substituting these expressions for and 0 D into ~f y 



'tf _ 

} M 



But gilberts = 0.1i7rNI ampere-turns 



so, 



NI = 2gBp _ 

OohrrMo 



2(0.1ii5)[il70 = 766 ampere-turns 

0 . Ml. 259) 



g = 0„lli5 cm. 

(measured) 

/*•* = 1.2^9 maxwell/NI-cm. 
Bp = hl70 gauss 
(measured) 



Bp= hl70 gauss was determined by use of a General Electric fluxmeter 
Model 32C2li8 which was calibrated at the National Bureau of Standards by the 
author for this purpose with the assistance of Mr. I. L. Cooter and A. R. 
Lindberg of the staff there on 2U August, 1951. Substituting in the equation 
for Kg and letting 

x ~ 0, Kgmin = 9.01 volt sec/inch 
x =0.032, Kgmax =12.6 volt sec/inch 



3.126 Electromechanical Coupling Constant 

Kjj, the electromechanical coupling constant, may be defined as 

follows: 

_ F (force in lbs) 

1 (differential current in amps) 

But F (0p -t— 0 a )^ from section 3.111 upon making the proper 

288S 



substitutions in the formula for force. 



( 0 p*b 0 a) 

Therefore, ${ = 288s where 0p = 13.9 kilomaxwells (measured) 

i 

0 a = 60 56 kilomaxwells (Fig. ^ ) 

S = gap area, 0.262 sq. in. 
i = 60 milliamperes 

Upon substituting the numerics and solving, is found to be: 

(13. 94-6. 56) 2 _ 92.8 lbs/ amp 

( 288) (0. 262) (0.060 ) 

3.2 PLOTS OF THEORETICAL IMPEDANCE AND PHASE AND AMPLITUDE CURVES 
3.21 Impedance Curves 

In section 3.11ii it was shown that the Laplace transform of the 
impedance was: 

2 - /g a - l[s 3 -b (f/M f R/L)s g + (k/M b fR/LM ■+- K b K m /LM)s b kR/LM ] 

1 s^ -b f s/M -b k/M 

or multiplying through by L, 

Z = Ls 3 ± (fL/M R)s ? b (kL/M b fR/M b KbKm/M)s +■ kR/M 
s 2 -b f s/M -t- k/M 

Making use of the properties of the Laplace transform, since 3 = 0 " -pjco 
with^'=0 J s = jco. Now, writing the impedance as a function of frequency 
by substituting jw for s, one has: 

Z(j») =/ Ae 1 (ow) _ L(» 3 ■+- (fL/M -+- R)(;w) 2 -b (kL/M -b fR/Mb KgK M /M)>tkR/M 
1 (»^bf(jw)/M -Hk/M 

or separating the real and imaginary parts 

Z(» = [R(k/M- M ? )-fI^ ? A]-b jw [L(k/M-^) b fR/M b ^/m] 

(k/M-« 2 ) -b j(f/M)w 

R in this case equals R c only. 

From this equation for impedance plots have been made in Fig. 13 of Z as a 
function of frequency. Also the variation of Z as the viscous damping term 



00 / 




/OP 2-00 300 4°° 

FIG. 13 IMPEDANCE CURVES FOR SEVERAL VALUES OF t 



TABLE II 



CALCULATION Diili* FOR TiiEJRETICiiL I-*U'ED>iNC& 



CURV 







f/ii 


= 0 


f/rt = 


47.5 


f/Ji = 


300 






y 


= o ! 


._J.= 


0.013 


..... J_= 


0.113 


f 


60 


L 


0 


Z 


0 


L 


0 


cps 


r^o/sec 

i 


otua^ decrees ! 


ohms 


decrees 


ohjis 


decrees 


0 


0 


1300 


0.0 


1300 


0.0 


1600 


0.0 


£0 


146 


5900 


53.0 






3050 


53.9 


40 


451 


5150 


69.4 










60 


376 


7450 


76.0 j 






7570 


75.1 


60 


50^ 


9970 


79.6 










100 


646 


15700 


31.9 : 


15700 


65.3 


15700 


79.9 


150 


754 


15900 


33.5 i 










140 


330 


10400 


64.7 






19400 


79.3 


160 


1003 


44500 


65.7 










160 


Ij.30 


33300 


36.9 ! 


3^700 


33.3 


59700 


70.5 


100 


1x01 


45600 


37 . 6 






35700 


63. 5 


533 


1457 


69600 


3o. 5 - 


63800 


7 3.6 


34500 


53.1 


510 


1340 


CXD 


-9 0.0 


148000 


7.7 


31600 


33.7 


540 


1330 


33500 


-37.3 


37300 


-53.0 


55500 


31.x 


^30 


1444 


7600 


-76.3 


7650 


-73.6 


19100 


19.6 


<~40 


1507 


5640 


(3 vO • ^ 


5760 


+ 61.7 


-u6100 


45.3 


444 


1 4 


5340 


+74 . 0 










^50 


1570 


6650 


76.0 










c60 


1630 


15300 


61.6 






16500 


65.4 


460 


1760 


17500 


34.1 










300 


io35 


41100 


35 . 1 


21300 


34.0 


25100 


33.1 


340 


4140 


47000 


3 6 . iv j 










360 


4x,60 


53400 


36.5 










400 


^510 


v^ooOO 


66.9 


33600 


36.0 


34500 


36.0 


450 


4 650 


39430 


67 . 4 ; 











is varied. Calculation data for these theoretical curves may be found in 
Table II. 



3.22 Freouency Response Curves 

In section 3.HU it was shown that the Laplace transform of the 
displacement per unit volt was: 
x _ Km/LM 

( f/M + R/L ) s 2 +- ( k/K -I- f R/LM h K b K m /LM ) s +- kR/LM 

or writing this as a function of frequency as explained in 3.21 

x( fa ) = K m /LH 

/«®g [kR/LM - (f/M -(- R/L)^ 2 ] -fj<o[k/M 4- fR/LM +- KgK M /LM - 

Plots have been made of x/tte g versus frequency for various values 

of Jo In Fig. lii curves are plotted considering the tube as a zero impedance 

source, i.e. r =0. Hie extremely large phase lags prohibit the use of 
P 

such a source for driving the motor in the servo loop outlined in Fig. 1. 

In Fig. 15 curves are plotted considering the plate resistance, r , to be 
equal to 60,000 ohms and it is obvious from a comparison of the change in 
the phase and amplitude characteristics why pentodes with their high plate 
resistances are used to drive the motor. Calculation data for these curves 
may be found in Table III. 

Plots have also been made of the displacement per unit ampere, 
x/i, versus frequency assuming a constant current source in Fig. 16. 
Calculation data may be found in Table IV. 



3.23 Transient Response 

To secure the transient response characteristic of the motor a 
step input voltage, K? is applied. That is. 



/xe g = f(t) — 



0, t<0 

k£ t>o 



j~r^7 




/OO ZOO 300 - 400 

FIG. 15 PLOTS OF DISPLACEMENT PER VOLT VERSUS FREQUENCE, R = 61,800 UHMS 




TABLE III CALCULATION DATA BOH PLOTS OF DISPLACEMENT PER UNIT VOLT AS h FUNCTION 

OF FREQUENCY 
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FIG. 16 PLOTS OF DISPLACEMENT PER AMPERE VERSUS FREQUENCY 




TABLE IV CALCULATION DATA FOR PLOT OF OUTPUT DIS LaCHHENT 
PER AJPERE VERSUb FREQUENCY 



f 

0 

20 
100 
140 
I 160 
bOO 
bio 

b20 

b30 

b40 

b60 

300 

400 



J =0.013 



x/i 

0 . 437 
0.442 
0.567 
0.793 
1.65 
4.45 
lb. 10 
4.13 
2.19 
1.39 
0.336 
0.413 
0.167 



0 

0 

0 

-1.3 

2.5 
6.7 
20. 5 
90.0 
153.9 
163.6 

174.5 
175.1 
177.*. 

178.5 



II 

o 

• 


113 


x/i 

- — - -f - 


' J» 


0.437 | 


0 


0.442 


0 


0.559 


-3 


0.764 


15.4 


l.~3 


36.4 


1.35 


67.0 


1 . 92 


90.0 


1.70 


lib. 3 


1.38 


IbS. c 


1.03 


140.1 ; 


0.739 


151.7 ! 


0.400 


16b. 7 ; 


0.165 


170.6 ' 



i 



Taking the Laplace transform, jue^ = K/s 
Substituting this for /ce in section 3.11U, we have 

/ o 

X Km jffs 

LMs 3 +- (fL -h MR)s 2 -H (kL + f R -h K B %)s +- kR 

or x «M fd 

s [_LMs 3 -h (fL -j-MR)s 2 -t- (kL +■ fR + K B % ) s +R ] \ j 

Dividing the numerator and denominator by kR, we have: 

X _ *mA R J K*) 

s RLM/kR)s 3 -I- (fL/kR 1 - M/k)s 2 -j- (L/R -t-r/k+- K B K M /kR)s +-l| ( J 

and substituting in the numerical values (values listed in appendix 

Table XVI under motor B) 

3T__ 9.6h x 10~ 6 l— 

s( 2 . 29 xlO _ 10 s 3 Hr 76.75xlO _3 s 2 -}- 6.07xl0 _ ^s l) 1 J 



Factoring the cubic in the denominator, 

9.61 x 10" 6 



X — 



J » 

s +- 2.95 x 10 3 J | (s + 203) +- jll90][(s-J-203) - J119ojl 



If now the inverse Laplace transform is taken, 

s-6 



x = 9.6hxl0" 



.-6 



1.000 - 0.l62e" 29 ^ 0t -t- 1.000e" 203t sin(ll90t -123.1° )J 



At t = 0, x =9.6UxlO (1.000 - 0.162 - 0.838)K 



and at t=oo, x = 9.6UxlO _6 K* 

If K=l, that is a unit step voltage is applied, 

x = 9.6hxlO“ 6 -1.56xlO- 29 ^ ot +- 9.6UxlO- 6 e~ 203t sin(1190t-123.1°) 



A plot of the transient response to a step input of one volt is shown in 
Fig. 17 and the calculation data in Table V. 




FIG. 17 TRANSIENT RESPONSE TO STEP INPUT 



TaLLE V CiiLCULa'IIOiJ DaTa I OR THEORETICAL TRaaEIaaT 
R^oPOadL TO a UMIT STEP IaPUT 

X=[9.64xl0" 6 -1.56xl0" 6 e" fc95Jt +J.64xl0“ 6 e" i: ' C) ' 3t siii( ll^Ot-lLS. T;] 
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U.O EXPERIMENT 

In this section the actual experimental techniques used to 
determine the constants derived analytically in section 3.0 on theory 
will be described. It will be noted that two motors, called A and B, 
were used. When the analytical section was begun it had been planned 
to use Motor A throughout both the theoretical and the experimental 
phases but in the course of the determination of the phase and amplitude 
characteristics the laminations composing the armature of Motor A separated 
and therefore a different motor had to be substituted. This new motor has 
been called Motor B and differs from Motor A only in that it is a separate 
physical entity subject to the usual variations in manufacturing processes. 

For the comparison of theoretical and experimental results a new set of 
theoretical curves for Motor B was computed. The values of Motor A have 
been left in the report for the purposes of comparison and to retain the 
analytical studies already completed. 

Iwl EXPERIMENTAL DETERMINATION OF ALL CONSTANTS 
h.ll Effective Mass 

The moment of inertia of a body can be determined experimentally 
from the following formula (see Ref. U, p. 72). 

where I_ = moment of inertia of 
a 9 

armature in gm-cm . 

T = period of oscillation sec/ cycle. 

2w = suspension width in cm. 

M a =mass of armature in grams. 

1 = length of suspension in cm. 

0 = acceleration of gravity, cm/sec^. 




